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PACS 12 . 20 . -m - Quantum electrodynamics 
PACS 42 . 50 . Ct - Ligth interaction with matter 
PACS 44 . 40 . +a - Radiative transfer in heat transfer 

Abstract. - We study the Casimir-Lifshitz force and the radiative heat transfer occurring be- 
tween two arbitrary bodies, each one held at a given temperature, surrounded by environmental 
radiation at a third temperature. The system, in stationary configuration out of thermal equi- 
librium, is characterized by a force and a heat transfer depending on the three temperatures, 
and explicitly expressed in terms of the scattering operators of each body. We find a closed-form 
analytic expression valid for bodies of any geometry and dielectric properties. As an example, 
the force between two parallel slabs of finite thickness is calculated, showing the importance of 
the environmental temperature as well as the occurrence of a repulsive interaction. An analytic 
expression is also provided for the force acting on an atom in front of a slab. Our predictions can 
be relevant for experimental and technological purposes. 
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Introduction. — The quantum and thermal fluctua- 
tions of the electromagnetic field result in a force between 
any couple of polarizable bodies. This dispersion effect, 
originally predicted by Casimir for two parallel perfectly 
conducting plates at zero temperature, and later gener- 
alized by Lifshitz and coworkers to real dielectric materi- 
als at finite temperatures [l], becomes relevant for bodies 
separated by less than few microns. Today, the Casimir- 
Lifshitz force plays a major role in all fundamental and 
technological issues occurring at these separations, becom- 
ing widely interesting, from biological systems to micro- 
electromechanics f2l|3l. 



The force shows two components, related to the purely 
quantum and thermal field fluctuations, respectively. The 
former dominate at short distances, and have been mea- 
sured in several configurations since early times. The lat- 
ter are relevant at large separations, with a weak total 
effect measured only recently [4j[5]. The Casimir-Lifshitz 
force, largely studied for systems at thermal equilibrium, 
has been considered also for systems out of thermal equi- 
librium [6U7], as done for the atom-atom, plane-plane and 
atom-surface force in [8 10 . More recently, other partic- 
ular configurations have been studied out of thermal equi- 
librium, involving infinitely thick planar and non planar 



bodies, and atoms 11 -15. Indeed, non-equilibrium ef- 
fects received a renewed interest since the recent discovery 
that systems driven out of thermal equilibrium may show 
new qualitative and quantitative behaviors, namely the 
possibility of a repulsive force, and of a strong force tun- 
ability 
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Such peculiar characteristics allowed a non- 
equilibrium system to be used for the first measurement 
of thermal effects, by exploiting trapped ultracold atomic 
gases (5]. 

Such a force, related to the correlations of the electro- 
magnetic field, shares a common formalism with radiative 
heat transfer [17], recently object of several investigation 
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for both fundamental and technological issues 

Due to a less direct formalism, calculations of forces and 
heat transfer for systems out of equilibrium have been per- 
formed only for some ideal configurations (atoms at zero 
temperature, infinitely extended surfaces, . . . ), and a gen- 
eral theory able to take into account arbitrary materials 
and geometries at different temperatures in an arbitrary 
thermal environment, is still missing. In this paper we 
derive such a theory, and we obtain a closed-form explicit 
expression for the most general problem of the force and 
heat transfer between two arbitrary bodies at two differ- 
ent temperatures placed in a thermal environment having 
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Fig. 1: (color online) Scheme of the system. The body 1 (2) is 
held at temperature T\ (T2), while the environment has tem- 
perature T3. 



a third (in general different) temperature. In addressing 
this problem we describe the two bodies by means of their 
scattering operators. This approach, successfully used to 
calculate the Casimir-Lifshitz force both at 
out 
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13 of thermal equilibrium, even if formally equiva- 



lent to a Green function technique, presents the advan- 
tage of requiring only single-body operators rather than 
the complete solution for the composite system. 

The physical system. — The system is made of two 
bodies, labeled with indexes 1 and 2, which we assume 
separated by at least an infinite plane. This assump- 
tion, introduced to exclude two concatenated bodies, is 
not necessary in approaches based on scattering theory 



20 , nonetheless it makes the formalism much easier to 
follow and it is verified practically in all the common ex- 
perimental configurations (e.g. plane-plane, sphere-plane, 
atom-surface). We are going to calculate the component 
of the forces acting on the two bodies along the axis per- 
pendicular to such a plane, to which we refer as the z axis, 
as well as the heat flux on each of them. A scheme of the 
system is shown in figure [T] where three distinct regions 
A, B and C are defined. The body 1 (2) is assumed to 
be at a local thermal equilibrium at temperature T\ (T 2 ), 
and on the two bodies impinges a thermal radiation at the 
third temperature T3. One can imagine a system where 
the two bodies are inside, and far from the surfaces, of a 
much bigger cell (the reservoir), held locally at thermal 
equilibrium at temperature T 3 . The three temperatures 
remain constant in time, so that the system assumes a 
stationary regime. If at least one of the two bodies is mi- 
croscopic, the validity of his hypothesis requires a more 
careful analysis: generally speaking it remains acceptable 
for time scales much shorter that the typical time of evo- 
lution of the microscopic body. 

The force F acting on any of the two bodies and the 
heat absorbed by it per unit of time can be evaluated 
calculating the fluxes 



(T(R,t)) Bym -dS, H 



(S(R, i)) S ym 



dS 

(1) 



of the symmetrized averages (AB) sym = (AB + BA)/2 
of the Maxwell stress tensor T and Poynting vector S in 
vacuum 



Tij (R, t) = e [EiEj + fBiBj - - (e 2 + c 2 B 2 



S(R,t) = e c 2 E x B 



(2) 



= x,y,z), through an arbitrary closed and oriented 
surface S enclosing the body. Choosing £ as a paral- 
lelepiped having two of its faces orthogonal to the z axis 
(and on two opposite sides of a given body) and letting 
the surface of these sides tend to infinity, one finds that 
the to component (to = x, y, z) F m of the force (the heat 
transfer) has non negligible contribution only from the dif- 
ference of the fluxes of the (T mz ) component of the tensor 
(the (S z ) component of the Poynting vector) on the two 
sides of the considered body. For the electric and magnetic 
fields appearing in eq. ^ we use the following mode de- 
composition 



E(R, t) = 2Rc 



+ OO 



2V 



E(R,w) (3) 



where 



E(R,c) = £ / ^e«*-^(k >W )^(k ) 



«,). (4) 



In this re presenta tion, known as angular spectrum repre- 
sentation 21 22 , a mode has amplitude £c(k, w), cor- 
responding to the frequency w, the transverse wavevec- 
tor k = (k x ,k v ), the transverse polarization p taking the 
values 1 (TE) and 2 (TM), and the direction of propa- 
gation along the z axis cf> = ±1 (with shorthand nota- 
tion (j> = ± in the polarization vectors and field ampli- 
tudes). In this approach, which proves to be convenient 
in our planar-like geometry, k z = \Jlo 2 jc 2 — k 2 is a de- 
pendent variable, and the three-dimensional wavevector 
is noted as = (k, (f>k z ). We have also introduced the 
polarization unit vectors, defined as e^ E (k, uS) = z x k 
and e^ M (k,w) = e^ E (k, uS) x K^ 1 where z = (0,0, 1) and 
A = A/ A. The analogous expression for the magnetic 
field can be directly deduced from Maxwell's equations. 
The field E^* in fig. [j]is the part of the field emitted by 
the body 7 (for 7 = 1,2), of the environmental field (for 
7 = 3), or of the total field in each region (for 7 = A, B, C) 
propagating in direction (f>. The fluxes of (T mz ) and (S z ) 
through a given plane z = z can now be explicitly ex- 
pressed in the following way 
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where 7 denotes the region in which z is located. Each 
of these contributions do not depend on the coordinate 
of the plane chosen to calculate the flux, provided that 
this surface is located inside a given region A, B or C. 
These quantities are expressed as a function of the matrix 
(jh)<f><f>' 



= 2n5{u - w')(p, k|C (7) ^' \p',k') 

(7) I 
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Fig. 2: Schemes for scattering operators. 

the two bodies, one recovers the force acting on the body 



1 along the z axis 19 20 



being the correlating function between two amplitudes of 
the total field propagating in directions <p an d <t>' , and 
associated to a couple of modes (k,p) and (k',p'). We 
have explicitly inserted the conservation of frequency, di- 
rect consequence of the time invariance (moving bodies are 
excluded from our calculation) characterizing our system. 
In eq. Q the correlators are written as matrix elements 
of a matrix C** : in our notation, the matrices are defined 
on the space (p, k) being p = 1, 2 and k£R 2 and thus the 
product of two matrices A and B is given by 



(p, k\AB\p', k') 



E 

p" 



2\JI 



d 2 k 



(27T) 

x (p",k"\B\p',k' 



2 {p,k\A\p",k" 



(8) 



At equilibrium. — For a system at thermal equilib- 
rium, the correlators are directly given by the fluctuation- 
dissipation theorem 



{E i (R,Lj)E*(R' : u')) sym = 2Tr5(LJ-u')-N(tj,T) 

u) (9) 

x ImGij(R, R',w), 

where the purely quantum and the thermal fluctua- 
tions are clearly distinguishable in the factor N(lj, T) = 
fiw[l/2 + n(u),T)}, with n{u,T) = (exp[fiw/(fc s T)] - l)" 1 , 
and Gij are the components of the Green function G as- 
sociated to the system of the two bodies, solution of the 
differential equation 



2 2 
^e( W ,R)jG(R,R>) = ^ 



Vr x V R - =Ve(w,R) G(R,R» = ^IS(R-R') 

(10) 

being I the identity dyad and e(w, R) the dielectric func- 
tion of the system. By expressing the Green function in 
terms of the reflection scattering matrices operators 1Z of 



F^m = -2ReTr(— N(u,T) 



(11) 



where = (1 - ftW+ftt 2 )-)- 1 and = (1 - 

72.( 2 ) _ 72-^ 1 ' )+ ) _1 . The definitions of the trace and the reflec- 



tion operators in ( 11 ) will become explicit in the following. 



Out of equilibrium. — For a system out of thermal 
equilibrium, the theorem (|9| is not valid, and the expres- 
sion of the correlators is not explicit in general. Neverthe- 
less, in the particular case of stationary non-equilibrium, 
this is possible by tracing back the knowledge of the cor- 
relators to the description of the fields emitted by each 
body alone and by the environment: this will be done by 
defining the operators describing the scattering produced 
by the presence of each body. To this end, taking into 
account a single body (say object 1), we consider an in- 
coming field coming from the left side (region A in figure 
[TJ), as shown in figure [2ja): this field produces, in general, 
a field on both sides of the body. 

We call the field produced in region A reflected and 
the one in region B transmitted, and we introduce the 
operators 1Z~ and T + relating the amplitudes of these 
outgoing fields to the amplitudes of the incoming field 
E^ n)+ (k,uj). Gathering in a vectorial notation all the 
incident, reflected and transmitted modes (for any p, k 
and u) E^+, E^- and E^ respectively, the defini- 
tions of the scattering operators read £?( re )~ = TZ~E^ + 
and i?( tr )+ = T + E^ m ' + . An analogous procedure defines 
1Z + and T~ (see figure |2fb)). If a body is at rest, the 
scattering process conserves the frequency and thus the 
matrix element (p, k, v\S\p', k! , uj 1 ) of any scattering oper- 
ator S is proportional to 2ttS(uj — u/). For convenience, 
we will work from now on with scattering operators S{uS) 
at a given frequency oj and defined in the subspace (p, k) . 
Before going further, it is useful to introduce a modified 
transmission operator which gives only the scattered 
part of the transmitted field. It is defined by the relation 
= l+T^ and, differently from the ordinary T^, it goes 
to zero in the limit of absence of the body. 
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To build the correlators, we need the expression of the 
total field in each region 7 = A, B, C of figure [lj which 
originates from the fields emitted by the bodies and the 
environment. We would like to express the total field in 
each region through the scattering operators. This can be 
simply obtained by summing up all the possible multiple- 
scattering events: for example the fields propagating in 
the two directions in region B are the solutions of 



E (B)+ = E (l)+ + j-(l) + £(3)+ + ft(l) + E (B)- 
E (B)- = E (2)- +7 -(2)- £ ;(3)- + 7^(2)-£;(B) + 



(12) 



The solutions of these equations as well as the analogous 
relations for regions A and C are straightforward and will 
not be given here explicitly. 

We are then left with the calculation of the correlators 
of the fields produced by the bodies and the environment. 
As far as the environment is concerned, it corresponds to 
a free bosonic field at temperature T 3 having correlators 



(if >nk, W )4? )0 V^'))s ym = S^-^N(u,T) 

Z€qC 



x Re(— J 5 pp ,(2Tr) 3 5(uj- u')5(k-k') 



(13) 



As for the bodies 1 and 2, our assumption that a local tem- 
perature can be defined for each one, and remains constant 
in time, reasonably leads to assume that the part of the 
total field emitted by each body is the same it would be 
if the body was at thermal equilibrium with the environ- 
ment at its own temperature. This hypothesis, already 
used in [5^ [TT[|T^[T6pT] , implies that the correlators of the 
field emitted by each body can still be obtained using the 
fluctuation-dissipation theorem (J9J) at its local tempera- 
ture, where the Green function is now associated to each 
body in absence of the other one. 

This procedure, together with the explicit connection 
between Green function and scattering operators (directly 
obtained from their definitions |23] ) allows us to obtain an 
explicit expression of the correlators of the field emitted by 
each body i = 1,2: in particular, for modes propagating 
in the same direction (0 = <j)') we obtain 



(£?W*(k J o;)£7«* t (k / jW / )) 



sym 



2e c- 



In these expressions, which represent a crucial intermedi- 
ate result, we have introduced the notation p(P w ' ew > = 
fc m n ( P w/ow) ( thig definition will be also used for m ^ 1 in 
the following), where n^ pw ' and n(° w ' are the projectors 
on the propagative and evanescent sector, respectively. 

The knowledge of the correlators (13), ( [l4| and (15) 
allows the evaluation of the flux of (T zz ) and (S z ) in each 
region. Finally, the differences of such fluxes ([l]) provides a 
closed-form analytic expression of the z component of the 
force and the heat transfer relative to body 1 expressed in 
terms of the three temperatures Ti, T 2 , and T3, and of the 
scattering operators of body 1 and 2: 

F(T U T 2 ,T 3 ) = Ficq) ( T i) + F(cq) (T2) +A2{Ti ^ T3) 



H(T 1 ,T 2 ,T 3 )=A 1 (T 1 ,T 2 ,T 3 ). 



(16) 



In eq. ( 16 ) we present the final result for the force as 



a sum of two contributions, the first term being a ther- 
mal average of the force F^ oq ^ (T) at thermal equilibrium 
given by eq. (11), at the temperatures of the two bodies 
, 12 . The two terms Ai and A 2 defined in 
( 16 1 can be collected as The force and heat trans- 



Xi and T 2 
eq 



fer on object 2 can be obtained from ( 16 ) by changing 



the sign, and by interchanging indexes 1 and 2, as well 
as + and — in its explicit expression. The term (17) is 
purely a non-equilibrium contribution, obeying the equal- 
ity A m (T, T, T) = 0. We remark that eq. §T7§ contains 
terms proportional to the transmission operators T^ 1 ^ 
and T^ 2 ' - , resulting from taking into account the finite- 
ness of objects 1 and 2, which were absent in previous in- 
vestigations concerning infinitely thick bodies 11-13. It 
is worth stressing that eq. 



(17) provides a finite value of 



the force and the heat transfer for any finite body 1. This 
property is not evident at first sight since the expression 
(17) contains also divergent terms. For instance, this is 



the case of the first term in the second line, proportional 
to the operator V^Z[- Since this operator is by defini- 
tion diagonal in the (k,p) basis, its trace is proportional 
to (2ir) 2 5(0), thus divergent. Moreover, we observe that 
this term diverges for any choice of the body 1, since it 
is formally independent on its scattering operators. Nev- 



ertheless, the quantity (17) remains finite due to peculiar 



(p, k| ^-p (p 1 w) - -Ji( i )4"p(P™)'ji(i)4> 1 i _ f(i)<l>'p(v™)f{i)<t>] 



while for (f> ^ <// we obtain 

(4^k )W )4^Vv)> s: 



2e n c< 



(p, k| f-H^p^ji*^ - 7-W^p(P w )-ft 



cancellations of several divergent terms. In order to show 
that, we must pay attention to all the terms which, in 
analogy with V^_{ , do not go to zero in absence of body 
_ 1. This is the case of the operators U^ 12 \ U {2 ^ and T<W 
\ 1 1 (0 = +, — ) ; having all the identity operator as a limit- 
ing value in absence of the body 1. By making use of 
the definition of the operator and of the properties 
N(uj,Ti)2Tr5{uj-Lj') C/( 12 ) = 1 + TZ^+TZ^-U^ (and analogous for U^) it 
is straightforward to show that all the terms not going to 
zero in absence of body 1 perfectly cancel each other, and 
only contributions proportional to either TZ^ 1 ^ or "J^ 1 )^ 
operators remain, leading indeed to a finite value of both 
(15) |the force and the heat transfer on any body 1. 
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A m (T 1 ,T 2 ,T 3 ) = (-l) m+1 h^{u; 2 - m [ r ^(A m (n^-,n^+) - (-l)"M m (7e«+, 7^ 2 >-)) + n^pWR^lT^RW" 

+(-i) m « m3 ^rl + ^3i[(-i) m s m (^ (1)+ ,^( 2 )-,r (1)+ ) - (^(pw)^(i)- 7 ,(pw) r (i)+t c/ (i2)t^(2)-t r (i)-t +hx .y 
+ [ni2 (v^ + n^-v^ - p^ftW-t - ^(2)-p(pw)^(2)-t^ + ^gr^-pi^r^^t (17) 

+n 13 ^ 2 )-T« + Pi7V« + ^ 

Equations (11) and (17) contain a trace, defined by the relation TtA(uj) = J2 P I (f^p Jo +0 ° ^ (Pj k|^4(o-') k) , the 
function riij — n(uj,Ti) — n(w, Tj) = 1,2,3), and the two supplementary functions 

A m (K {1)+ ,K {2) -) = U (12) (V ( ^ ] - K^+V^K^ +ft (1) V c 1 w) - 7>i7^ (1)+t ) 

x u {12) Uvlr ] + {-iTn^v^n^- +n {2) -^v£ w) + {~i) m v^ w) n (2) -), 

(18) 

B m (n^+ \n^- ,r (1)+ ) = [/( 12 )r (1)+ 7> p TT (1)+t c/ (12)t 

x (><p w) + {-i) m nW-^v^nW- +n (2) -^v£ w) + (-i) m v^ w) n {2) - 



Some applications. — We will now apply eq. (17) 



to two specific systems, namely atom-surface and slab- 
slab configurations. As anticipated before, although the 



derivation of the main result ( 17) benefits from the choice 
of the plane-wave basis, it can be applied to any cou- 
ple of bodies (separated by a plane) even not satisfy- 
ing translational invariance. An interesting candidate to 
check the validity of this feature is the atom-plane sys- 
tem. In this case, the scattering operators associted to 
the atom can be deduced by describing the atom as an 
induced dipole d(uj) — a(w)E(R J 4, uj) proportional to the 
uj component of the electric field calculated at the atomic 
position R^i = (r^, za) through the atomic dynamical po- 
larizability a(uj) (isotropy for the atomic polarizability has 
been assumed) . Writing the field produced by the induced 
dipole and projecting it on the plane wave basis as done 
in 24 we obtain the following expression for the reflection 



and transmission atomic operators (for 



<k, P KH|kV; 



iuj 2 a{uj) 



2enc 2 fc z 



(^(k, W )-e/(k>)) 



x exp[i(k' - k) • r^] exp[-i0(fc 2 + k' z )zA\ 
iuj 2 a(uj) 



(k,p|T A »|k>'; 



2enc 2 fc, 



■(e£(k jW ).e*(k / , 



x exp[i(k' — k) ■ r^] exp[— i(f>(k z — k' z )zA\- 



(19) 



These expressions clearly show that the wavevector k and 
the polarization p are not conserved in the atomic scat- 
tering process, as a consequence of the absence of transla- 
tional invariance. These operators can be inserted in eqs. 



(11) and (17) in order to obtain the equilibrium and non- 
equilibrium contributions to the force: to this aim, only 
the terms up to the first order in 1Z and T are kept, co- 
herently with a second-order perturbative approach with 



respect to the electric charge (first order with respect to 
atomic polarizability). As shown for example in 24 , this 
gives back the known result for the force at thermal equi- 
librium. As far as the non-equilibrium contribution is con- 
cerned, taking the atom on the left (za < 0) of the surface 
z = 0, a simple calculation yields 



A 2 (7i, T 2 ,T 3 



rim 



(Y.I 



+ 00 



4ir 2 e c z ^ Jo 



dojuj 2 a(us) 



"23 



dkk(\ Pp \ 2 + \T p \ 2 -l) 



+ I ' dkk(e+ ■e-)(n 32Pp e- 2ik ' ZA + n 13 p* p e 2lk * ZA ) 



f+oo 

nyi I dk k(e 



£-\ r ,*„-2ik z z A 



e P )P P e 



(20) 



where p p and t p are respectively the reflection and trans- 
mission Fresnel coefficients associated to a planar slab and 
the dependence of all the quantities inside the integral 
on co and k is implicit. We remark the fact that in eq. 



(20), the polarizability a(oj) is for an atom at tempera- 
ture T\ . The first term in the square bracket in eq. ( 20 ) is 



a distance-independent contribution already discussed in 
[l0p6 . The second and the third terms are both distance- 
dependent: the former depends on the propagative sector, 
whereas for the latter only evanescent waves contribute. 
We have verified that, taking a groun-state atom (Xi = 0) 
and assuming that T 2 and T 3 are such that atomic exci- 
tation can be excluded (which amounts to replace a(u>) 
by its static value cc(0)), we recover the result of 16], ob- 
tained using a different approach. 

We consider now a second example for which we show 
quantitatively the new physical features produced by these 
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Fig. 3: (color online) Pressure acting on a 2/im thick slab par- 
allel to a lOOO^im thick slab (see text). Lines: equilibrium pres- 
sures at T = OK (solid), 300 K (dashed), 600 K (dash-dotted). 
Symbols: non-equilibrium pressures, T3 = OK (circles), 300 K 
(diamonds), 600 K (crosses), with Ti = 300 K and T 2 = OK in 
(a)-(b) and Ti = T 2 = 300 K in (c)-(d). 



terms. We calculate the non-equilibrium pressure on a 
2 fim thick slab (body 1, fused silica) parallel to a 1000 /im 
thick slab (body 2, silicon). In figure [3]ja)-(b) we show 
the case T x = 300 K and T 2 = OK, whereas in[3];c)-(d) 
T\ = T 2 = 300 K: in both cases T3 takes the three values 
K, 300 K and 600 K. The figure shows that in both cases 
(a)-(b) and (c)-(d) the non-equilibrium pressure can sig- 
nificantly differ from the equilibrium counterpart at any 
of the three temperatures involved. Moreover, both for 
equal and unequal Ti and T 2 , the choice T 3 = pro- 
duces a repulsive force starting around 6 /Ltm of distance 
between the plates. This is particularly remarkable in the 
case Ti = T 2 = 300 K, showing that the environmental 
temperature may play an important role for objects of fi- 
nite thickness, qualitatively modifying the behavior of the 
force. 

Conclusions. — We have derived a general expres- 
sion for the Casimir-Lifshitz force and for the radiative 
heat transfer for systems out of thermal equilibrium, valid 
for bodies having arbitrary shape and dielectric function. 
Depending on the bodies and on the environmental tem- 
peratures, the force and heat transfer present several in- 
teresting degree of freedom. Due to its generality, our 
results allow a straightforward study of the force and heat 
transfer for systems involving bodies whose scattering ma- 
trices are analytically known (atoms, cylinders, spheres 
and also the investigation of most 



and slabs) p23 25 26 



general bodies by a numerical evaluation of the scatter- 
ing matrix. In particular, the heat transfer expression will 
allow to obtain more accurate estimations useful for past 
and future experiments as well as for technological appli- 



cations such as solar cells. The force, which we calculated 
explicitly for an atom in front of a slab and numerically for 
two parallel slabs, can be significantly affected by thermal 
non-equilibrium with the environment, presenting transi- 
tion from attractive to repulsive behaviors at distances of 
few microns. 
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